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Abstract 

Background: In equine dentistry periodontal diseases, especially periapical inflammation, are frequently occurring 
problems. Anachoresis is believed to be the most common cause for the development of such disorders. 
Nevertheless, there is still no substantiated explanation why settlement of pathogen microorganisms occurs in 
equine periodontal tissues. It is expected that excessive strains and stresses occurring in the periodontal ligament 
(PDL) during the horse's chewing cycle might be a predisposing factor. In this study this assumption was examined 
by finite element (FE) analyses on virtual 3-D models of equine maxillary and mandibular cheek teeth, established 
on the basis of ljCT datasets. Calculations were conducted both under conditions of closing and power stroke. 

Results: Results showed a uniform distribution of low stresses and strain energy density (SED) during closing stroke, 
whereas during power stroke an occurrence of high stresses and SED could be observed in the PDL near the 
alveolar crest and in periapical regions. 

Conclusion: The concentration offerees during power stroke in these specific areas of the PDL may cause local 
tissue necrosis and inflammation and thus establish a suitable environment for the settlement of microorganisms. 

Keywords: Finite element analysis, Horse, Periodontal ligament, Tooth, Chewing cycle, Periodontal disease, 
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Background 

Dentition plays a very important role in equine digestion. 
Horses do not have a forestomach system like ruminants, 
which enables a microbiological breakdown of plant cell 
walls before reaching the resorptive small intestine [1]. 
Thus, cheek teeth are the only tools for releasing nutri- 
tive cell contents by grinding the herbivore diet [1]. An 
effective disruption of forage is therefore necessary for a 
sufficient energy supply [2] . 

One of the main problems in equine oral health are 
painful periodontal diseases [3-5], hampering normal 
masticatory action [6]. Especially in geriatric horses the 
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occurrence of periodontal disorders is exceptionally high 
with described prevalence of about 60% [6]. 

Periodontal diseases can be subdivided into those pro- 
ceeding from the gingival sulcus into the periodontal 
space and those occurring exclusively in the periapical 
region. 

Recent studies have shown that periapical infections 
can be caused either by periodontal spread, infundibular 
caries or occlusal fissure fractures [7-9]. Nonetheless, re- 
markably the aetiology of most cases of the periapical in- 
fection (68%) remained unexplained [7]. For these cases 
anachoresis is suggested to be causative [7-9]. Anachore- 
tic infections depend on a suitable environment for 
microbiological settlement in terms of necrotic tissue 
areas [10]. 

We propose that excessive strains and stresses occur- 
ring in the periodontal ligament (PDL) during the horses 
chewing cycle could cause pathological tissue changes 
and thus be a predisposing factor for equine periodontal 
diseases. 

Ltd. This is an Open Access article distributed under the terms of the Creative 
:ommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and 
iginal work is properly cited. 
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In this study, computer-based finite element (FE) cal- 
culations of stresses and strain energy densities (SED) 
occurring in the PDL of different- aged teeth during load 
bearing phases of the equine masticatory cycle (closing 
stroke and power stroke) were conducted. The results 
will hopefully give useful contributions for understanding 
the etiology of periodontal diseases. These calculations 
were generated on the basis of our previously established 
finite element models of equine maxillary and mandibu- 
lar cheek teeth [11,12]. 

Methods 

3-D models 

The establishment of virtual 3-D models of different 
aged mandibular and maxillary equine cheek teeth (tooth 
age group A: 0-5 years; age group B: 6-15 years; age 
group C: > 15 years), suitable for FE-simulations, was 
conducted on the basis of (iCT-datasets [12]. In all 3-D 
models the tooth was assumed to consist of one 
homogenous material, while when constructing the jaw 
bone a separation into cortical alveolus (lamina dura), 
cancellous bone and compact bone was made. The PDL, 
connecting the tooth with the surrounding cortical al- 
veolus, was built uniformly around the tooth except for 
periodontal parts around the apical region of the tooth 
and under the root-bi-/trifurcation, which were segmen- 
ted as different, modified materials [11]. Surface and vol- 
ume meshes, necessary for finite element analyses, were 
generated out of the 3-D models as described by Lupke 
et al. (2010) [12]. The quality of all meshes was tested 
afterwards by a quality measuring function of the com- 
puter program COMSOL Multiphysics (version 3.4, 
COMSOL AB, Stockholm, Sweden) to ensure that mesh 
topology did not negatively affect the FE solutions. 

Finite element analyses 

Finite element simulations of the horses chewing cycle 
were performed for three defined age groups (A-C) using 
the computer program COMSOL Multiphysics (version 
3.4, COMSOL AB, Stockholm, Sweden). 

For calculating stresses and SED defined values of 
Youngs modulus and Poissons ratio (values describing 
the elastic behavior of materials) had to be allocated to 
each structure in the 3-D meshes (Table 1). Youngs 
moduli and Poissons ratios for the tooth [13-15] and 
for the bony structures [16,17] as well as the Poissons 
ratio for the PDL [18,19] were shown to be uniform in 
several studies and were therefore taken from litera- 
ture. In this study both the tooth and all bony struc- 
tures were assumed to have linear elastic material 
properties. 

The published values for Youngs modulus of the 
PDL vary in a wide range dependent on species and 
measuring methods [20-22]. Thus this value was 



Table 1 Values of Poisson's ratio and Young's modulus 
used in FE calculations 





Poisson's 
ratio 


Young's modulus 
(MPa) 


References 


Tooth 


0.3 


20000 


[13-15] 


Compact bone 


0.3 


20000 


[16,17] 


Cancelous bone 


0.3 


1000 


[16,17] 


PDL 


0.45 


3.1 (age group A) 


[11] 






2.9 (age group B) 








2.6 (age group C) 





specifically determined by intrusive movement experi- 
ments for each of the used models as described else- 
where [11]. Due to the viscoelastic properties of the 
PDL [23-25] a non-linear elastic behavior was proven. 
Nevertheless, for later FE calculations, especially 
regarding the quality but not quantity of results, an al- 
location of a linear but age-dependent Youngs modu- 
lus to the PDL was carried out. 

All structures in the model were assumed to be 
isotropic. 

The amount of masticatory forces which were applied 
in the FE simulations was measured and calculated in 
previous studies [26,27]. For simulations of the closing 
stroke a force of 350N [26] was applied perpendicular to 
a plane described by the alveolar crest (Figure 1). This 
plane was assumed to be perpendicular to the tooths 
longitudinal axis. 

The power stroke was simulated using a force of 
1250N. The 3-D orientation of the force vector, applied 
in FE calculations, was obtained from analyses of the 
equine chewing cycle [1,28]. According to these data 
the masticatory forces during the power stroke cause 
an axial intrusion of the tooth into its alveolar socket 
both in maxillary and mandibular teeth, but also an 
important sideward movement (shifting) of the tooth 
as a consequence of the apparent friction. Due to the 
movement of the mandible from lateral to medial, 
maxillary teeth undergo shifting into a palatal direction 
while mandibular teeth undergo shifting in a buccal 
direction. This phenomenon was simulated in FE cal- 
culations by applying the force at an angle of 45° to 
the plane through the alveolar crest. In maxillary teeth 
this force was orientated in a palatal direction while in 
mandibular teeth it was orientated in a buccal direc- 
tion (Figure 1). Efforts were made to ensure that the 
force was applied to the tooth along the transverse 
ridges of the occlusal surface, which determine the 
pathway of the mandible during occlusal contact of op- 
posing teeth [1]. 

A fixation of the outer buccal and palatal/lingual sur- 
faces of the compact bone was chosen as a necessary 
physical requirement in all calculations. 
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(1) Closing stroke 



(2) Power stroke 



350 N 



1250 N 




Figure 1 Force application for closing stroke and power stroke. Forces were applied perpendicular to a plane through the alveolar crest for 
closing stroke (la) and at an angle of 45° to this plane for power stroke (2a) . Force application occurs along the transverse ridges for power 
stroke (b). 



Results 

Depiction of finite element simulations 

Stresses, occurring in the periodontal ligament in occluso- 
apical direction (s z normal stress in z-direction, i.e. direc- 
tion of intrusion) were calculated both for closing and 
power stroke. Stresses occurring in linguo-buccal or 
bucco-palatal directions (s x normal stress in x-direction, i. 
e. shifting direction) were additionally computed for power 
stroke because this direction was also assumed to be 
placed under considerable force during this phase of the 
chewing cycle. Stresses in a mesio-distal direction (s y nor- 
mal stress in y-direction) were not depicted because the 
amount of these stresses was insignificant compared to the 
s x and s z stresses, and irrelevant regarding loading during 
power stroke. Furthermore, analyses of the SED were con- 
ducted for closing and power stroke. 

Due to the shifting of teeth during power stroke, for 
every model a distinction should be made between the side 
to which the tooth is pressed (the shifting side, i.e. palatal 
side of maxillary teeth, buccal side of mandibular teeth) 
and the opposite side (the shifting- opposed side, i.e. buccal 
side of maxillary and lingual side of mandibular teeth). To 
be consistent this distinction is used in the depiction of 
results both for the power stroke and closing stroke (see 
Figures 2, 3, 4, 5 and 6). 



Finite element simulations 

The simulation of the closing stroke showed a uniform 
distribution of predominantly tensile stresses in an axial 
direction (s z stresses) within the PDL for the teeth of age 
group A (young) and B (middle) (Figure 2). In age group 
C (old) stress distribution appeared irregular, showing 



tensile but also compressive stresses especially in the 
mandibular model. 

Although the dimensions of these calculated s z stresses 
during closing stroke increased with tooth age, values 
were far below those calculated for the power stroke (see 
below). 

The distribution of strain energy densities within the 
PDL during closing stroke (Figure 3) was quite uniform 
for all three age groups except for the results of the man- 
dibular tooth of the old age group, where accumulations 
of SEDs were calculated analogously to the computed 
stress accumulations in Figure 2. An increase of SEDs 
with increasing tooth age was also noticed. 

During power stroke a specific distribution pattern of 
compressive and tensile stresses was obtained. In maxil- 
lary cheek teeth tensile stresses both in an occluso-apical 
direction (s z stresses) and in linguo-buccal or bucco-pal- 
atal directions (s x stresses) occurred in the PDL in 
regions near the alveolar crest on the buccal side (shifting- 
opposed side) and in apical regions of the palatal side 
(shifting side). In these models with a linear-elastic 
assumed PDL, compressive stresses (s z , s x ) were calculated 
in regions near the alveolar crest on the palatal side and in 
apical regions of the buccal side (Figures 4, 5). 

Calculations on mandibular cheek teeth showed similar 
stress distributions: Tensile stresses occurred near the al- 
veolar crest on the lingual side (shifting-opposed side) 
and in apical regions of the buccal side (shifting side), 
while compressive stresses were calculated near the al- 
veolar crest on the buccal and in apical regions on the 
lingual side (Figures 4, 5). Stresses occurring during 
power stroke were higher than those occurring during 
closing stroke. 
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Figure 2 Stresses occurring in the PDL in occluso-apical direction (sz normal stresses) during closing stroke. A load of 350N was applied 

along the tooth's longitudinal axis. Diagram shows FE results for the three different age groups (A-C) on both the shifting side (1) and side 

where shifting side is averted (2). Tensile stresses are demonstrated by positive scale values (yellow to red), compressive stresses by negative scale 

values (green to blue). Dimension of stresses increases with increasing tooth age (please note the different age-related scales), 
v ) 



The highest amplitudes of SED during power stroke 
occurred in the PDL of both maxillary and mandibular 
cheek teeth especially in regions near the alveolar crest 
and also in periapical regions (Figure 6). 

Thus, for all age-groups, in both maxillary and man- 
dibular cheek teeth, major stress or SED locations could 
be observed near the alveolar crest and in periapical 



regions during power stroke. Furthermore, an increase in 
stresses and SED with increasing dental age could be 
noticed. 

The fact, that on each side of the tooth (shifting side 
and shifting-opposed side) both compressive and tensile 
stresses occurred, showed that the tooth was not uni- 
formly shifted in one direction during the power stroke 
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AT B 




sz normal 
stress 
(MPa) 

Max 4.5 





B2 



C1 



C2 



sz normal 
stress 
(MPa) 

Max 6 






C1 



C2 



Figure 4 Stresses occurring in the PDL in occluso-apical direction (s z normal stresses) during power stroke. A load of 1250N was applied 
at an angle of 45° to a plane through the alveolar crest. Diagram shows FE results for the three different age-groups (A-C) on both the shifting 
side (1) and side where shifting side is averted (2). Tensile stresses are demonstrated by positive scale values (yellow to red), compressive stresses 
by negative scale values (green to blue). Dimension of stresses increases with increasing tooth age (please note the different age-related scales). 



but that it is rotated in its alveolus around a center of ro- 
tation. This could be determined as being approximately 
in the middle of the reserve crown. 

Discussion 

The horses masticatory cycle is subdivided into three 
different phases: opening stroke, closing stroke and 



power stroke [1,28]. Contact of maxillary and mandibu- 
lar cheek teeth and thus forces in the PDL occur during 
the closing and power stroke. Therefore finite element 
calculations were performed for each of theses phases. 
During the closing stroke a mean force of 350N [26] 
occurs in occluso-apical direction along the longitudinal 
axis of the tooth, which simply causes an intrusion of the 




sx normal 
stress 
(MPa) 

Max 3 




sx normal 
stress 
(MPa) 

Max 4.5 



C1 



C2 



(MPa) 
Max 6 







C1 



C2 



Figure 5 Stresses occurring in the PDL in bucco-palatal / linguo-buccal direction (s x normal stresses) during power stroke. A load of 
1250N was applied at an angle of 45° to a plane through the alveolar crest. Diagram shows FE results for the three different age-groups (A-C) on 
both the shifting side (1) and side where shifting side is averted (2). Tensile stresses were demonstrated by positive scale values (yellow to red), 
compressive stresses by negative scale values (green to blue). Dimension of stresses increases with increasing tooth age (please note the different 
age-related scales). 
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tooth into its alveolus. In our calculations this event was 
represented by a force application perpendicular to the 
alveolar crest. 

The power stroke causes a grinding of the lower jaw 
teeth upon the upper jaw teeth with a mean load of 
1250N [26]. This load is divided into an occluso-apical 
orientated vector, representing the compression of the 
maxillary and mandibular teeth against each other, and a 
linguo-buccal (mandibular teeth) or bucco-palatal (max- 
illary teeth) orientated vector, describing the shearing ac- 
tion and friction between the teeth due to the sideward 
movement of the mandible. The ratio of the force of fric- 
tion between two bodies and the force pressing them to- 
gether is described by the coefficient of friction (COF). 

As no experimental data about the COF during the 
equine power stroke were available, the assumption of 
both vectors having the same absolute values (COF = 1) 
was made to primarily describe primarily the quality of 
the power stroke. The shifting of the tooth in its alveolus 
should be simulated but the occurring axial intrusion 
should not be neglected either. Thus, the force of 1250N 
was applied to the alveolar crest at an angle of 45°. 

In this study, both stresses and strain energy density 
were calculated. The strain energy density was computed 
due to the knowledge that this elastic strain energy, 
stored in the material per volume unit, is one parameter 
which is related to cell reaction and remodeling pro- 
cesses in bones [29]. Since other studies also describe 
cell reactions of tendons (for example, in terms of a re- 
lease of PGE 2 ) due to repetitive mechanical loading 
[30,31] and explain cell movement resulting from strain 



energy density as a mechanical stimulus [31,32], the 
strain energy density is also expected to be correlated to 
tissue responses within the periodontal ligament [33]. 
The threshold SED, which results in cell reaction, could 
be variable [29]. Nonetheless, it is expected that the 
higher the SED the higher the probability of cellular re- 
action [33,34]. 

In using the preliminary determined Youngs moduli 
for the PDL of different age-groups, disregarding the 
orthotropic behavior of the PDL (different material para- 
meters in different orthogonal directions) must always 
be taken into consideration when interpreting the calcu- 
lated stresses and SED. The determined Youngs moduli 
were actually only precisely determined for dental move- 
ment/intrusion along the tooths longitudinal axis but 
not for displacements in other directions (buccolingual/ 
mesiodistal). Because of the age-dependent histology of 
the periodontal ligament, in terms of changing density 
and orientation of collagen fibers [35,36], the PDLs 
Youngs Modulus is expected to be different for other 
directions than the axial one. Since stress calculation 
depends on the Youngs modulus, the absolute amount 
of stresses occurring in buccopalatal/linguobuccal direc- 
tion (s x stresses) should be considered with caution. 

Further, the allocation of a linear Youngs modulus (in- 
stead of a non-linear Youngs modulus) to the PDL with 
its viscoelastic properties instead represents quite a sim- 
plification. It was proved while applying a linear Youngs 
modulus to the PDL in finite element analyses that cal- 
culated compressive forces were higher than those oc- 
curring in vivo, whereas tensile forces were definitely 
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underestimated by the model [13]. This is why predic- 
tions concerning the absolute quantity of stress dimen- 
sions and the occurrence of computed compressive 
stresses should be interpreted with caution. Nevertheless, 
such simplifications would not influence the general con- 
clusions of these simulations concerning major localiza- 
tions of occurring stresses and strain energy density. 
Thus, it can be expected that the finite element calcula- 
tions of the closing and power stroke in this study deliv- 
ered useful and reliable results in view of the 
understanding of the etiology of periodontal diseases and 
principles of dental treatment. 

During closing stroke the simulations show a quite 
uniform distribution of low stresses and SED compared 
to the power stroke. In age-groups A (0-5 years) and B 
(6-15 years) mainly tensile stresses were calculated, this 
meeting the expectations of a conduction of chewing 
forces predominately by tension of collagen fiber bundles 
[13,23]. Computed compressive stresses in the old age 
group (>15 years) were calculated because of difficulties 
concerning an exact force application in axial direction. 
Due to the anatomical curvatures in cheek teeth, their 
longitudinal axes are difficult to determine. For this rea- 
son, and in terms of reproducibility, the longitudinal axis 
in these simulations was assumed to be perpendicular to 
a plane described by the alveolar crest. Such inevitable 
inaccuracy in axial force application might lead to an in- 
advertent shifting of the tooth in a distinct direction. 
Such changes in force application lead to greater displa- 
cements in older teeth due to their smaller anchorage 
surface in the alveolar socket. This explains the com- 
puted compressive stresses and SEDs in the mandibular 
model of the old teeth. Basically, the closing stroke is 
probably not responsible for the development of peri- 
odontal disorders. 

The most important phase considering the etiology of 
periodontal diseases might be the power stroke, where 
quite high stresses and SED occur in regions near the al- 
veolar crest and around the root tips. 

Forces around the alveolar crest might contribute to 
the development of gingivitis by periodontal damage 
resulting in gingival pockets, which would deliver a per- 
fect entrance for microorganisms. It is known that peri- 
odontal damage in terms of fiber rupture occurs under 
stresses above 3.8N/mm 2 [37]. In FE calculations this 
value was exceeded during power stroke especially in the 
PDL of the old age group. Therefore, on the basis of our 
results it can be expected that especially old horses are 
affected by gingivitis, a fact that is already described in 
literature [5]. The highest stress and SED concentrations 
in the PDL of the old age group can be explained by the 
steadily decreasing anchorage surface for the periodontal 
fibers in aging and shortening teeth. The PDL is the 
most important tissue for the attenuation and 



conduction of chewing forces [38]. Thus, it became clear 
that the exposure of the small PDL of old teeth under 
chewing forces is higher than that of the large PDL of 
young teeth. 

Further, occurring forces might directly contribute to the 
development of diastemata which are inevitably found to be 
associated with periodontal disease [39]. The development 
of diastemata results from tooth movements. Such move- 
ments are proved to be initiated by deformations of the 
PDL subsequently stimulating the process of alveolar bone 
remodeling [40]. Thus, it can be expected for the equine 
dentition that occurring SED in the PDL initiates biomech- 
anical and cellular reactions leading to bone remodeling 
resulting in diastemata. In our study it was shown that SED 
in the PDL increases with dental age. This leads to the as- 
sumption that the prevalence of periodontal diseases should 
increase with age, too. This, indeed, coincides with other 
studies describing periodontal diseases especially in older 
horses [5,39]. 

The concentration of forces in apical regions is known to 
cause local trauma and tissue necrosis [41], which are suit- 
able conditions for the settlement of microorganisms [10]. 
High concentrations of stresses and SED occurring in the 
old age group could explain the appearance of periapical 
infections in old horses. According to our results the 
described high incidence of periapical infections in young 
horses [39] cannot be explained by stresses and SED result- 
ing from chewing forces. In young horses the periapical 
regions are still in an early odontological phase. They are 
surrounded by a dental organ and dental sack. In these 
areas proliferative processes inevitably cause massive bone 
remodeling resulting in the well-known phenomenon of 
eruption cysts. It is assumed that the osteoclastic action, 
causing eruption cysts, is initiated by forces within the sur- 
rounding tissues. These forces are presumably generated by 
periodontal tissues but not by masticatory forces [42]. Thus, 
the hypothesis that high stresses in the periapical region 
promote the development of anachoresis is supported. 

In this study one mean force of 1250N was taken for all 
calculations in all different age groups to improve reprodu- 
cibility of single results. However, studies on masticatory 
forces in horses [26,43] show that forces acting on single 
teeth vary within a wide range (875N in Triadan position 06 
to 1970N in Triadan position 11). The actual forces are 
influenced by two factors: First, by the position of the teeth, 
with forces increasing in caudal direction, and second, by 
the Curve of Spee, which is the curve which connects the 
cusps of the maxillary teeth and is tangent to the mandibu- 
lar condyle. Due to the flattening of the Curve of Spee with 
increasing age, forces acting on single teeth decrease with 
age [43]. The mean force of 1250N used in our study 
reflects the forces in Triadan position 09 (for all age groups) 
quite well. Forces acting in Triadan postions 10 and 11 are 
known to be higher for all age groups [26]. Thus, in using a 
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mean force of 1250N in our study, calculated stresses and 
SED might be overestimated, presumably in premolars of 
old horses, or - more importantly - even underestimated, 
presumably for Triadan positions 10 and 11. 

Remarkably, this assumption is mirrored by the described 
frequent location of diastemata and following periodontal 
disorders between second and third molars [44]. 

FE calculations with defined forces to match individual 
tooth positions should be performed in future studies to 
describe the clinical difference of periapical infections 
subject to Triadan-position. 

Conclusion 

Occurring forces in the equine PDL are highest during the 
power stroke. We propose that high stresses, strains and 
SED near the alveolar crest and in periapical regions predis- 
pose to the development of periodontal diseases and peria- 
pical infections by causing local trauma and thus 
establishing an optimal environment for microbiological 
settlement. From FE calculations, showing high forces oc- 
curring in the PDL of old teeth, a higher prevalence of peri- 
odontal diseases in older horses can be concluded. 
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